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qq | This paper proposes a new method for estimating the contribution from different 

CN ' ionospheric regions to the response of total electron content variations to the solar flare, 

based on data from the international network of two-frequency multichannel receivers 
\ of the navigation GPS system. The method uses the effect of partial shadowing of 

the atmosphere by the terrestrial globe. The study of the solar flare influence on the 



Oh- 



atmosphere uses GPS stations located near the boundary of the shadow on the ground 
kjpj, in the nightside hemisphere. The beams between the satellite-borne transmitter and the 

^ | receiver on the ground for these stations pass partially through the atmosphere lying 

O ■ in the region of total shadow and partially through the illuminated atmosphere. The 

analysis of the ionospheric effect of a powerful solar flare of class X5.7/3B that was 
recorded on July 14, 2000 (10:24 UT, N22W07) in quiet geomagnetic conditions (Dst = 
O -10 nT) has shown that about 20% of the TEC increase correspond to the ionospheric 

region lying below 100 km, about 5% refer to the ionospheric E-region (100-140 km), 
about 30% correspond to the ionospheric Fl-region (140-200 km), and about 30% to 
' regions lying above 300 km. 
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1 Introduction 

The enhancement of X-ray and ultraviolet (UV) emission that is observed during chromospheric flares on 
the Sun immediately causes an increase in electron density in the ionosphere. These density variations are 
different for different altitudes and are called Sudden Ionospheric Disturbances, SID (Davies,1990; Donnelly, 
1969). SIDs are generally recorded as the short wave fadeout, SWF (Stonehocker, 1970), sudden phase 
anomaly, SPA (Ohshio, 1971), sudden frequency deviation, SFD (Donnelly, 1971; Liu et al., 1996), sudden 
cosmic noise absorption, SCNA (Deshpande and Mitra, 1972), sudden enhancement/decrease of atmospher- 
ics, SES (Sao et al., 1970). Much research is devoted to SID studies, among them a number of thorough 
reviews (Mitra, 1974; Davies, 1990). 

Highly informative technique is the method of Incoherent Scatter (IS). The Millstone Hill IS facility 
recorded a powerful flare on August 7, 1972 (Mendillo and Evans, 1974a). The measurements were made in 
the height range from 125 to 1200 km. The increase of local electron density N e made up 100% at 125 km 
altitude and 60% at 200 km. 

Using the IS method Thome and Wagner (1971) obtained important evidence of the height distribution 
of the increase in N e at the time of the May 21 and 23, 1967 flares. A significant increase of N e was recorded 
in the E-region, up to 200%, which gradually decreases in the F-region with the increasing height, down to 
10-30%, and remains distinguishable up to 300 km. The earliest increase of N e begins in the E-region, and 
at higher altitudes it is observed with a delay which is particularly pronounced at F-region heights. 

A sudden increase in total electron content (TEC) can be measured using continuously operating radio 
beacons installed on geostationary satellites. On August 7, 1972, Mendillo et al. (1974b) were the first to 
make an attempt to carry out global observations of the solar flare using 17 stations in North America, 
Europe, and Africa. The observations covered an area, the boundaries of which were separated by 70° in 
latitude and by 10 hours in local time. For different stations, the absolute value of the TEC increase Al 



varies from 1.8 • 10 16 to 8.6 • 10 16 el • m~ 2 , which corresponds to 15-30% of the TEC. Investigations revealed 
a latitudinal dependence of the TEC increase value. At low latitudes, it was higher compared with high 
latitudes. Besides, the authors point out the absence of a connection between the TEC increase value and 
the solar zenith angle. 

The advent and evolution of a Global Positioning System (GPS) and also the creation on its basis of 
widely branched networks of GPS stations (at least 900 sites at the August of 2001, the data from which are 
placed on the Internet) opened up a new era in remote ionospheric sensing. High-precision measurements 
of the TEC along the line-of-sight (LOS) between the receiver on the ground and transmitters on the GPS 
system satellites covering the reception zone are made using two-frequency multichannel receivers of the 
GPS system at almost any point of the globe and at any time simultaneously at two coherently coupled 
frequencies /i = 1575.42 MHz and / 2 = 1227.60 MHz. 

The sensitivity of phase measurements in the GPS system is sufficient for detecting irregularities with an 
amplitude of up to 10 3 — 10 4 of the diurnal TEC variation. This makes it possible to formulate the problem 
of detecting ionospheric disturbances from different sources of artificial and natural origins. The TEC unit 
(TECU) which is equal to 10 16 el • m -2 and is commonly accepted in the literature, will be used throughout 
the text. 

Afraimovich (2000a); Afraimovich et al. (2000b, 2001a,b) developed a novel technology of a global 
detection of ionospheric effects from solar flares and presented data from first GPS measurements of global 
response of the ionosphere to powerful impulsive flares of July 29, 1999, and December 28, 1999. Authors 
found that fluctuations of TEC are coherent for all stations on the dayside of the Earth. The time profile of 
TEC responses is similar to the time behavior of hard X-ray emission variations during flares in the energy 
range 25-35 keV if the relaxation time of electron density disturbances in the ionosphere of order 50-100 s is 
introduced. No such effect on the nightside of the Earth has been detected yet. 

Afraimovich ct al. (2001c) and Leonovich et al. (2001) suggested a technique for estimating the iono- 
spheric response to weak solar flares (of X-ray class C). They obtained a dependence of the ionospheric TEC 
increase amplitude (during the solar flare) on the flare location on the Sun (on the central meridian distance, 
CMD). For flares lying nearer to the disk center (CMD < 40°), an empirical dependence of the ionospheric 
TEC increase amplitude on the peak power of solar flares in the X-ray range was obtained (using data from 
the geostationary GOES-10 satellite). 

This paper is a logical continuation of the series of our publications (Afraimovich, 2000a; Afraimovich et 
al., 2000b, 2001a, b,c, Leonovich et al., 2001) devoted to the study of ionospheric effects of solar flares, based 
on data from the international GPS network. 

A limitation of the GPS method is that its results have an integral character, as a consequence of which 
from measurements at a single site it is impossible to determine which ionospheric region makes the main 
contribution to the TEC variation. The objective of this study is to develop a method which would help 
overcome (at least partially) this problem. 

2 Method of determining the shadow altitude ho over 
the ground 

The method uses the effect of partial "shadowing" of the atmosphere by the terrestrial globe. Direct beams 
of solar ionizing radiation from the flare do not penetrate the region of the Earths total shadow. GPS stations 
located near the shadow boundary on the ground in the nightside hemisphere are used to investigate the 
solar flare influence on the ionosphere. The LOS for these stations pass partially through the atmosphere 
lying in the total shadow region, and partially through the illuminated atmosphere. The altitude over the 
ground at which the LOS intersects the boundary of the total shadow cone, will be referred to as the shadow 
altitude h - 

Fig. 1 schematically represents the formation of the cone of the Earths total shadow (not to scale) in the 
geocentric solar-ecliptic coordinate system (GSE): the axis Z is directed to a north perpendicularly planes 
of an ecliptic, the axis X - on the Sun, the axis Y is directed perpendiculary to these axes. For definition of 
the shadow altitude ho it is necessary to know coordinates of a cross point C of the LOS and the shadow 
boundary. 

The primary data are the geographical coordinates of station GPS on the Earth (Fig.l; a point P): 
an elevation angle and azimuth of LOS on a satellite GPS, toward the north clockwise, for the time (UT) 
corresponding to the phase of solar flare maximum in the X-ray range. These coordinates are converted to 



the Cartesian coordinate system where the Cartesian coordinates of the GPS station on the ground and the 
coordinates of the subionosphcric point (at 300 km altitude) arc calculated. Next, we use the geocentric 
solar-ecliptic coordinate system following the technique reported by Sergeev and Tsygancnko (1980). To 
determine the coordinates of the point C we solve a system of equations: the equation of cone (of total 
shadow), and the equation of a straight line (LOS) specified parametrically. After that, from the resulting 
point C we drop a perpendicular to the ground and calculate its length (Fig. 1, line ho). The value of ho, 
thus obtained, is just the shadow altitude. 



3 Method of determining the TEC increase in the 
ionosphere using data from the global GPS net work 

This paper exemplifies an analysis of the ionospheric effect of a powerful solar flare of class X5.7/3B recorded 
on July 14, 2000 (10:24 UT, N22W07) under quiet geomagnetic conditions (Dst = -10 nT). The time profile 
of soft X-ray emission in the range 1-8 A (GOES-10 data) at the time of the flare is presented in Fig. 2a. 

To determine the TEC increase in the ionosphere we used the data from the international GPS network. 
The GPS technology provides a means of estimating the TEC variations Io(t) on the basis of TEC phase 
measurements made with each of the spatially separated two-frequency GPS receivers using the formula 
(Calais and Minster, 1996): 

Io(t) = 40.308 ff-"fi i{LlXl ~ L " X2) + C ° nSt + nL]: (1) 

where LiAi and L2X2 are the increments of the radio signal phase path caused by the phase delay in 
the ionosphere (m); L\ and L 2 stand for the number of complete phase rotations, and Ai and A 2 are the 
wavelengths (m) for the frequencies /1 and / 2 , respectively; const is some unknown initial phase path (m); 
and nL is the error in determining the phase path (m). 

Input data used in the analysis include series of the oblique value of TEC Io (t) , as well as corresponding 
series of elevations 9 and azimuths of LOS to the satellite. These parameters are calculated using our 
developed CONVTEC program to convert standard (for the GPS system) RINEX-files received via the 
Internet. Input series of TEC Io(t) are converted to the vertical value following a well-known technique 
(Klobuchar, 1986). 
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where Re is Earth's radius; and h max is the height of the ionospheric F2-layer maximum. 

Variations of the regular ionosphere, and also trends introduced by the motion of the satellite are elim- 
inated using the procedure of removing the trend defined as a polynomial of the third order on a given 
temporary interval. 

Figs. 2b and 2d present the typical time dependencies of the vertical TEC I(t) for sites GPS WDLM 
(PRN02, shadow altitude ho = 17 km) and LEEP (PRN07, shadow altitude ho = 586 km). Time dependen- 
cies of the TEC AI(t) response, with the trend removed for these series, arc presented in Figs. 2c and 2e, 
respectively. 



4 Results and discussion 

The TEC response to the solar flare was analyzed for 45 GPS stations. Detailed information about the GPS 
stations and analysis results is summarized in Table 1: names of GPS receiving stations (Site), number of 
the GPS satellite from which the signal is received (PRN), shadow altitude above the ground (ho), absolute 
increase of TEC A J, relative increase of TEC (A/(i)/A/oo(t), and geographical coordinates of GPS stations 
(latitude, longitude). The increase of TEC A/00 corresponds to the amplitude of the TEC increase measured 
at the station lying at the shadow boundary on the ground (ho = 0). 

Fig. 3b illustrates examples of time dependencies of the TEC AIo response for LOS to the satellite which 
intersect the boundary of the shadow cone at different heights ho during the solar flare of July 14, 2000. 
Fig. 3b (left) presents the values of these altitudes, and (right) the names of corresponding stations. For a 



better visualization, the dependencies are drawn by lines of a different thickness. It should be noted that 
the response remains pronounced when the shadow altitude exceeds significantly the electron density peak 
height in the ionosphere. For station GUAM (PRN26, height of the shadow boundary ho — 885 km), the 
response amplitude exceeds the background oscillation amplitude by more than a factor of 2. 

It is evident from Fig. 3 that the wave phase (time of the response maximum) is different at different 
altitudes ho- On the one hand, this phenomenon can be caused by the interference of the response with back- 
ground fluctuations; on the other, this can be due to the fact that at different heights different wavelengths 
of ionizing radiation are observed, which, in turn, can have independent time characteristics. 

The dependence of the absolute TEC increase on the altitude ho for all the cases under consideration is 
plotted in Fig. 4a. The dependence of the relative TEC increase AI(t)/AIoo on the altitude ho during the 
solar flare is shown as a percentage in Fig. 4b. The TEC increase A7oo(i) corresponds to the amplitude of 
the TEC increase measured at the station lying at the shadow boundary on the ground (ho = 0). 

Fig. 4b suggests that about 20% of the TEC increase correspond to the ionospheric region lying below 
100 km, about 5% refer to the ionospheric E-region (100-140 km), about 30% correspond to the ionospheric 
Fl-region (140-200 km), and about 30% to regions lying above 300 km. We found that a rather significant 
contribution to the TEC increase is made by ionospheric regions lying above 300 km. 

The estimate obtained is consistent with the findings reported by Mendillo and Evans (1974a); Mcndillo 
et al. (1974b). The authors of the cited references, based on investigating the electron density profile in the 
height range from 125 km to 1200 km using the IS method, concluded that about 40% of the TEC increase 
during the powerful flare on August 7, 1972, correspond to ionospheric regions lying above 300 km. However, 
Thome and Wagner (1971), who used the IS method to investigate the ionospheric effects from two others 
powerful solar flares, pointed out that an increase in electron density associated with the solar flare was 
observable to 300 km altitude only. This difference can be explained by the fact that each particular solar 
flare is a unique event which is characterized by its own spectrum and dynamics in the flare process. 

Acknowledgements. Authors are grateful to V.G. Mikhalkovsky for his assistance in preparing the English 
version of the manuscript. This work was done with support from both the Russian foundation for Basic 
Research (grant 00-05-72026) and RFBR grant of leading scientific schools of the Russian Federation No. 
00-15-98509. 

References 

Afraimovich, E. L., GPS global detection of the ionospheric response to solar flares, Radio Sci., 35, 
1417-1424, 2000a. 

Afraimovich, E. L., Kosogorov, E. A., and L. A. Leonovich, The use of the international GPS network as 
the global detector (GLOBDET) simultaneously observing sudden ionospheric disturbances, Earth Planet. 
Space, 52, 1077-1082, 2000b. 

Afraimovich, E.L., Altyntsev, A.T., Kosogorov, E.A., Larina, N.S., and L. A. Leonovich, Detecting of 
the Ionospheric effects of the solar flares as deduced from global GPS network data, Geomagnetism and 
Aeronomy, 41, 208-214, 2001a. 

Afraimovich, E.L., Altyntsev, A.T., Kosogorov, E.A., Larina, N.S., and L. A. Leonovich, Ionospheric 
effects of the solar flares of September 23, 1998 and July 29, 1999 as deduced from global GPS network data, 
J. Atm. Solar-Terr. Phys., 63, 1841-1849, 2001b. 

Afraimovich, E.L., Altyntsev, A.T., Grechnev, V.V., and L. A. Leonovich, Ionospheric effects of the solar 

flares as deduced from global GPS network data, Adv. Space Res., 27, 1333-1338, 2001c. 

Calais, E., and J. B. Minster, GPS detection of ionospheric perturbations following a Space Shuttle 

ascent, Geophys. Res. Lett, 23, 1897-1900, 1996. 

Davies K., Ionospheric radio, Peter Peregrinus, London, 1990. 

Deshpande, S. D. and A. P. Mitra, Ionospheric effects of solar flares, IV, electron density profiles 
deduced from measurements of SCNA's and VLF phase and amplitude, J. Atmos. Terr. Phys., 34, 255- 
259, 1972. 

Donnelly, R. F., Contribution of X-ray and EUV bursts of solar flares to Sudden frequency deviations, J. 
Geophys. Res., 74, 1873 - 1877, 1969. 

Donnelly, R. F., Extreme ultraviolet flashes of solar flares observed via sudden frequency deviations: 
experimental results, Solar Phys., 20, 188-203, 1971. 



Klobuchar, J. A., Ionospheric time-delay algorithm for single-frequency GPS users, IEEE Transactions 
on Aerospace and Electronics System, AES , 23(3), 325-331, 1986. 

Leonovich, L.A., Altynsev, A.T., Afraimovich, E. L., and V. V. Grcchnev, Ionospheric effects of the solar 



flares as deduced from global GPS network data. LANL e-print archive, 2001, http: / / arXiv.org/ abs/physics/OllOOC . 
Liu, J. Y., Chiu, C.S., and C. H. Lin, The solar flare radiation responsible for sudden frequency deviation 
and geomagnetic fluctuation, J. Geophys. Res., 101, 10855-10862, 1996. 

Mendillo, M., and J.V. Evans, Incoherent scatter observations of the ionospheric response to a large 
solar flare, Radio Sci, 9, 197-203, 1974a. 

Mendillo, M., Klobuchar, J. A., Fritz, R. B., da Rosa, A.V., Kersley, L., Yeh, K. C, Flaherty, B. J., 
Rangaswamy, S., Schmid, P. E., Evans, J. V., Schodel, J. P., Matsoukas, D. A., Koster, J. R., Webster, A. R., 
and P. Chin, Behavior of the Ionospheric F Region During the Great Solar Flare of August 7, 1972, J. 
Geophys. Res., 79, 665-672, 1974b. 

Mitra, A. P., Ionospheric effects of solar flares, D.Reidel, Norwell, Mass., 1974. 
Ohshio, M., Negative sudden phase anomaly, Nature, 229, 239-244, 1971. 

Sao, K., Yamashita, M., Tanahashi, S., Jindoh, H., and K. Ohta, Sudden enhancements (SEA) and decreases 
(DSA) of atmospherics, J. Atmos. Terr. Phys., 32, 1567-1573, 1970. 

Sergeev, V. A., and N. A. Tsyganenko, The earth's magneto sphere. Results of researches on the 
international geophysical projects. 'Nauka', Moscow, (in Russian), 1980. 

Stonehocker, G. H., Advanced telecommunication forecasting technique in AGY, 5th., Ionospheric fore- 
casting, AGARD Conf. Proc, 29, 27-31, 1970. 

Thome, G.D, and L. S. Wagner, Electron density enhancements in the E and F regions of the ionosphere 
during solar flares, J. Geophys. Res., 76, 6883-6895, 1971. 



Fig. 1. Scheme for the determination of the shadow altitude ho over the ground in the GSE 
system. GPS - GPS navigation satellite; P - GPS station on the ground; C - intersection point of 
the transmitter - receiver LOS and the shadow boundary. 



Table 1. Parameters of response TEC during solar 
flare of July 14, 2000 



Site 


PRN 


ho 


AI 


AI : A/oo 


latitude 


longitude 


GPS 




km 


TECU 


% 


degree 


degree 


STB1 


2 





1.1 


100 


44.7 


272 


KAYT 


21 





1.06 


96.3 


13.9 


120 


STL4 


2 


15 


1.1 


100 


38.6 


270 


WDLM 


2 


16 


1.07 


97.2 


44.6 


264 


SLAI 


2 


21 


1.09 


99.0 


41.9 


266 


GUS2 


2 


65 


1.0 


90 


58.3 


225 


PRDS 


9 


80 


1.1 


100 


50.8 


245 


CORD 


1 


92 


1.0 


90.9 


-31.5 


295 


WILL 


2 


107 


0.9 


81.8 


52.2 


237 


YAR1 


17 


112 


1.0 


90.9 


-29.0 


115 


PLTC 


2 


140 


0.83 


75.4 


40.1 


255 


VCIO 


7 


141 


0.8 


72.7 


36.0 


260 


TMGO 


2 


147 


0.72 


65.4 


40.1 


254 


PATT 


7 


147 


0.71 


64.5 


31.7 


264 


DSRC 


2 


150 


0.7 


63.6 


39.9 


254 


LKWY 


9 


152 


0.71 


64.5 


44.5 


249 


PERT 


30 


155 


0.66 


60 


-31.8 


115 


NANO 


2 


175 


0.54 


49.0 


49.2 


235 


WHD1 


2 


185 


0.5 


45.4 


48.3 


237 


SEAW 


2 


194 


0.41 


37.2 


47.6 


237 


UCLU 


2 


195 


0.5 


45.4 


48.9 


234 


SEAT 


2 


195 


0.49 


44.5 


47.6 


237 


LIND 


9 


205 


0.46 


41.8 


47 


239 


SEAW 


9 


206 


0.42 


38.1 


47.6 


237 


GOBS 


2 


221 


0.44 


40 


45.8 


239 


SATS 


2 


222 


0.34 


30.9 


46.9 


236 


GWEN 


2 


227 


0.32 


29.0 


45.7 


238 


LMUT 


2 


233 


0.4 


36.3 


40.2 


248 


LKHU 


4 


234 


0.39 


35.4 


29.9 


264 


AZCN 


7 


280 


0.35 


31.8 


36.8 


252 


FTS1 


7 


280 


0.3 


27.2 


46.2 


236 


SHIN 


2 


339 


0.29 


26.3 


40.5 


239 


TUNG 


9 


344 


0.28 


25.4 


40 


241 


SMYC 


9 


412 


0.3 


27.2 


36.3 


244 


JAB1 


17 


413 


0.3 


27.2 


-12.6 


132 


UCLP 


2 


483 


0.31 


28.1 


34.0 


241 


CVHS 


9 


520 


0.27 


24.5 


34.0 


242 


LEEP 


7 


585 


0.28 


25.4 


34.1 


241 


PMHS 


7 


590 


0.25 


22.7 


33.9 


241 


UCLP 


7 


590 


0.26 


23.6 


34.0 


241 


CSDH 


7 


594 


0.25 


22.7 


33.8 


241 


YBHB 


26 


685 


0.24 


21.8 


41.7 


237 


RIOG 


13 


704 


0.25 


22.7 


-53.7 


292 


PTSG 


26 


719 


0.26 


23.6 


41.7 


235 


GUAM 


26 


885 


0.23 


20.9 


15.5 


154 
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Fig. 2. Time profile of soft X-ray emission in the range 1-8 A (GOES 10 data) during the solar flare 
of July 14, 2000 - (a). Time dependencies of TEC variations for the stations: WDLM (PRN02, 
shadow altitude h = 17 km) - (b), and LEEP (PRN07, shadow height h = 586 km) - (d). These 
same time dependencies upon subtracting the trend, respectively, - (c, e). 
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Fig. 3. Time profile of soft X-ray emission in the range 1-8 A (GOES 10 data) - (a). Examples of 
time dependencies of the TEC response AI(t) for LOS intersecting the boundary of the shadow 
cone at different altitudes /i during the solar flare of July 14, 2000 (for a better visualization, the 
dependencies are drawn by lines of different thicknesses); the column at the left shows the values 
of these heights, and the column at the right shows the names of the corresponding stations - (b). 
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Fig. 4. Dependence of the absolute TEC increase AI(t) on the shadow altitude ho during the solar 
flare of July 14, 2000; the dashed line on the plot corresponds to the root-mean-square deviation 
of background oscillations in the absence of a flare - (a). The dependence of the relative TEC 
increase AI(t) / A/oo {t) on the shadow altitude ho during the solar flare in percent - (b). The TEC 
increase A/oo (*) corresponds to the amplitude of the TEC increase measured at the station lying 
on the shadow boundary on the ground ho = 0. 



